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1, INTRODUCTION 

Compression consists of reducing the physical size of an information block (the number of bits) 
needed to represent the characteristics of an image (eliminating redundancies) to enable the reduction of storage 
cost and the rapid transmission of data [1]. The current methods of compression of still and moving images 
(video) are based on the wavelet transform; this latter had a huge success in the field of image processing such 
as image compression due to its ability to compress the image energy on a small number of coefficients 
allowing efficient coding of the image [2, 3]. 

Traditional multiresolution decompositions appear to be a restricted and limited category [4]. In order 
to overcome this problem, new transforms better adapted to the representation of the images have been 
proposed. In this work, a new method based on the quincunx wavelet transform coupled with the SPIHT-Z 
progressive encoder had introduced to compress different types of images such as medical, satellite, and natural 
images in order to prove the effectiveness of the proposed method for compression with low bit rates.. 


2. QUINCUNX TRANSFORM 

The dyadic separable analysis requires three families of wavelets, this is sometimes considered as a 
disadvantage, on the other hand the expansion factor between two successive scales is 4, which may seem high. 
It 1s possible to remedy these two problems, but on the cost of the loss of the separability of the filters and 
therefore of a slightly higher algorithmic complexity. An analysis has been particularly studied and found 
practical applications; it is the so-called "quincunx" analysis. Figure 1 illustrates this type 
of decomposition [5]. We will see that the expansion factor is only 2 between two successive resolutions and 
only one family of wavelets is necessary [6]. 


1 1 
In this case, the expansion matrix equals M “| : 
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The transformation of the mesh (lattice) is done according to the following diagram: 
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Figure 1. Quincunx wavelet decomposition 


This matrix generates a quincunx lattice in two dimensions. The column vectors of this matrix form a 
base of this lattice. The volume of the associated elementary cell is 2 Figure 2. This same lattice also comes 


from the matrix [7, 8] 
1 -l 
M'= 
1 1 





Quincunx grid Elementary cell 


Figure 2. Example of the quincunx grid and the elementary cell 


It can be seen that the dilation pitch is 0 on each direction and the geometry of the mesh obtained 
justifies the name given to this multiresolution analysis [9]. 
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As the quincunx sampling reduces the image size by a factor of 2, then the associated filter bank has 
two channels Figure 3. The low-pass filter H reduces the resolution by a factor of V2 . The coefficients of 


the wavelet correspond to the output of the high- filter G 
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Figure 3. The perfect reconstruction filter bank with the quincunx sampling 


Where H and G (respectively H andG ) are the transfer functions of the synthesis filters 
(respectively analysis). These filters correspond to the bases of decomposition and reconstruction of the 
wavelets [10]. The perfect reconstruction tells us that there is no loss of information during the iterations of 
analysis. In our case the analysis will be orthogonal, 1.e., the analysis and synthesis filters are the same. 


The underlying quincunx wavelet is given by 


ValX=J2 LD galiilb,(M X-ii) (1) 


neZ 
Where g, [7] is the spatial version of the filter G, la]. 


The advantage of this approach is that filters for small values of @ are almost isotropic (Figure.4), that is, the 
coefficients of the obtained wavelet transform have no preferential orientation, and this reflects the quality of 


\S 





% sae 
SN “a 
S i. 
ae a Wl 
iil 


Ula 





(a) a= /2 


Figure 4. The perfect The isotropy of the filter H, = (e! e for different values of the order parameter a 


3. THE SPIHT CODING ALGORITHM 
3.1. The classic SPIHT algorithm 

The SPIHT algorithm (Set Partitioning In Hierarchical Trees) [12] performs a recursive partitioning 
of the tree so as to determine the position of the significant coefficients in the progeny of the considered 
coefficient. It is based on the same concepts: progressive coding by bit plane and use of hierarchical 
dependencies between the coefficients of the different subbands. However, a new dependency protocol 
between the coefficients 1s defined. The coefficients of the lower frequency sub-band (LL3 in Figure 5) are 
grouped by four, and for each group, the offspring is as follows: One of the four coefficients (LL3 in Figure 5) 
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does not allow descendants, while the other three have four descendants each. For the other subbands, each 
coefficient has four descendants [13] 
Three sets of coordinates are used to present the new coding method : 
a. D(x,y): Set of all the descendants of the pixel (x,y). 
b. O(%,y) : Set of coordinates of the direct descendants of the pixel (x,y). 
c. L(x,y): Set defined as being D(x,y)-O(x, y). 

The SPIHT algorithm uses three lists of significant coefficients (LSP), insignificant coefficients (LIP), 
and insignificant sets (LIS). The entries for each of the lists are the (x, y) coordinates in the image matrix, 
which represent individual pixels in the case of LIP and LSP lists, and the set of D(x, y) or L(x, y) descendants 
in the case of the LIS list [14]. 





Figure 5. Repartition image used by SPIHT algorithm 


3.2. The modified SPIHT algorithm (SPIHT-Z) 

Figure 6 illustrates the tree of the SPIHT-Z algorithm according to the quincunx decomposition of the 
image. The representation of the spatial tree of the SPIHT-Z is as follow: The coefficient (x, y) of the subband 
(H6) points to two sons (a set of 2 x 1 adjacent pixels) (x-Lx,y+ Ly) and (x- Lx+1,y+Ly). The coefficients of 
the subband (H5) point to four descendants (x-Lx,y+ Ly) ,(x- Lx+1,y+Ly), (x+Lx,y- Ly+1) and (x+ Lx+l,y- 
Ly), of coefficients of the subband (H4) and so on until the last subband H1 as shown in Figure 6. 





Figure 6. Repartition image used by SPIHT-Z algorithm 
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The SPIHT-Z algorithm can be summarized as follow 

a. If the number i of the sub-band (H,) is an even integer, the pixel (x, y) points to two pixels of the sub- 
band (H;_,) having the coordinates (x — L,., y + Ly) and (x — L, + 1,y + L,,) as shown in Figure 7.(a). 

b. If the number of iterations i of the subband (H,) is an odd integer, the pixel (x, y) of the sub-band (H;) 
points to two pixels of the subband (H;_,) having the coordinates (x + L,,y — Ly +1) and (x +L, + 
1, y — L,,) as shown in Figure 7.(b). 


Lx 





(b) 


Figure 7. SPIHT-Z algorithm pattern 


This pattern which starts increasing from subband (H_1i) to higher subband (H_(-3)) draws a “Z” 
symbol which gives the appellation “SPIHT-Z” as shown in Figure 5. As the SPIHT algorithm, the SPIHT-Z 
considers three lists: LSP, LIP and LIS. 


4. QUALITY EVALUATION PARAMETERS 
The commonly parameter used in image compression is the mean square error MSE. This variable 


(error) is defined by the mean square between the pixel (i, J ) of the original image / (i J ), and the pixel (, J ) 


of the reconstructed image / (i 7 ) [15, 16]. 


MSE = awn j)- (i, i) (2) 


m=0 n=0 


The peak signal to noise ratio [17, 18] 
pies | 
PSNR =10.log,0 (tif [dB] (3) 
MSE 


We then evaluate a new paradigm for assessing the quality of medical images, the similarity index 
compares the brightness, contrast, and structure between each pair of vectors, where the index of structural 
similarity (SSIM) between two signals x and y is given by the following expression. 


SSIM (x, y) = (x, y)e(x, y)s(x, y) 


Finally, for the purpose, we need a single global measure of overall image quality that is given by 
the formula 


MSSIM(I, i)= = ssa I, (4) 
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M is the total number of local windows in the image. The values MSSIM exhibit greater consistency with the 
visual quality. [19, 20] 

The Visual Information fidelity parameter (VIF) quantifies Shannon information that is shared 
between the reference and distorted images with respect to the information contained in the reference image 
itself. VIF test 1s then evaluated as [21-23] 


_ yy 1(CI;F)|s4) 


VIF = ric 1(CJ;EJ|s/) 


(5) 


Where, [(X;Y|Z) is the conditional mutual information between X and Y, conditioned to Z; js isa 
realization of S j for a particlar image, the index j runs through all sub-bands in the decomposed image. 


5. RESULT 

In this paper, the quincunx wavelet transform (QWT) coupled with the SPIHT and SPIHT-Z coders 
had applied. a filter order parameter a=3 and an iteration number A=6 for wavelet quincunx decomposition. 
Tables 1, 2 and 3, given in appendises A, B and C, show respectively the compression results of 20 medical 
images, 20 satellite images and 20 natural images of size 512 x 512. All considered images are 255 gray-level. 
the bit rate varies from 0.1 to 0.5 bpp, and the image evaluation parameters PSNR, MSSIM and VIF are used 
to compare QWT CDF9/7 + SPIHT-Z, QWT CDF9/7 + SPIHT and DWT (with lifting structure [24-26]) 
CDF9/7 + SPIHT. 

From these results, it can be noted that the proposed algorithm QWT CDF9/7 + SPIHT-Z gives better 
image compression quality compared to the other algorithms (QWT CDF9/7 + SPIHT, DWT (Lifting) CDF9/7 
+ SPIHT), and this can be observed in the important values of PSNR = 46.83, MSSIM = 0.98, VIF = 0.74 for 
medical images; PSNR = 31.3, MSSIM = 0.81, VIF = 0.39 for satellite images; and PSNR = 37.77, MSSIM = 
0.92, VIF = 0.55 for natural images asb shown in Figure 11. 

In order to confirm these results, the averages of the evaluation parameters (PSNR, MSSIM, and VIF) 
of the 20 test images had calculated. Figures 8, 9 and 10 show the average curves of medical, satellite, and 
natural images. From these results, the proposed algorithm QWT CDF9 / 7 + SPIHT-Z gives good image 
compression quality especially at low bit rates. 
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Figure 8. Variation of average PSNR, MSSIM and VIF (over 20 Medical images) 
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Figure 9. Variation of average PSNR, MSSIM and VIF (over 20 satellite Images) 
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Figure 10. Variation of PSNR, MSSIM and VIF (over 20 natural Images) 
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PSNR=47.41dB, MSSIM=0.98 PSNR=31.3dB, MSSIM=0.81 PSNR=37.77dB, MSSIM=0.92, 


6. 


VIF=0.78 VIF=0.39 VIF=0.55 


Figure 11. Medical, satellite, and natural images compressed by CDF9/7(QWT) + SPIHT-Z 


CONCLUSION 
In this work, the efficiency of the quincunx wavelet transform QWT has been improved by introducing 


a modified SPIHT (SPIHT-Z) encoder that exploits the transition between all subbands of the quincunx wavelet 
transform (QWT). The obtained results are very satisfactory compared to other methods. The proposed 
algorithm is better adapted to low bit-rate image compression and especially in the medical field. 
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